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Entropy production, energy loss and currents in adiabatically rocked thermal 

ratchets 
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Abstract: We study the nature of currents, input energy and entropy production in different 
types of adiabatically rocked ratchets using the method of stochastic energetics. The currents 
exhibit a peak as a function of noise strength. We show that there is no underlying resonance or 
synchronisation phenomena in the dynamics of the particle with these current peaks. This follows 
from the analysis of energy loss in the medium. We also show that the maxima seen in current as 
well as the total entropy production are not directly correlated. 
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I. INTRODUCTION 



The subject of noise induced transport has attracted 
much theoretical as well as experimental interest for the 
past few years Q, 0, 0] • The motivation for such a 
study stems from the challenge to develop models to 
explain the reliable unidirectional transport observed in 
biological systems amidst a very noisy environment in 
the absence of bias. Systems that combine the asym- 
metry and nonequilibrium fluctuations to generate sys- 
tematic motion in the absence of a macroscopic bias 
are termed as ratchets or Brownian motors. In thermal 
equilibrium the principle of detailed balance prohibits 
any net particle current in the system. Hence a net cur- 
rent in the absence of any bias can appear only as a 
consequence of the interaction of the particle with its 
noisy nonequilibrium environment. Thus it is possible 
to extract energy from the random fluctuations and put 
it into use. These ratchet systems are information en- 
gines analogous to the Maxwell's demon which extract 
work out of bath at the expense of an overall increase 
in entropy (or entropy production) 0, [(J. There are 
several ways in which one can incorporate nonequilib- 
rium effects arising out of the irreversible interaction of 
the system with its external surroundings. This has led 
to various types of ratchets, namely, flashing ratchets, 
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rocking ratchets, frictional ratchets, time asymmetric 
ratchets, etc |2|. Extensive studies have been carried 
out on the nature of current and their possible rever- 
sals as a function of various physical parameters. These 
studies are found to be useful in identifying proper mod- 
els for biological motors and also to develop machines 
at the molecular scales including nanoparticle separa- 
tion devices Q. 

The subject of the energetics of Brownian motors or 
ratchets has developed into an entire subfield of its own 
right 0, 01 • A general framework has been developed 
wherein, the compatibility between the Langevin and 
the Fokker-Planck formalisms used for various types of 
ratchets or motor models and the laws of thermodynam- 
ics have been proved Q . Using this framework one can 
readily calculate various physical quantities like the ef- 
ficiency of energy transduction, energy dissipation (hys- 
terisis loss), entropy (entropy production), input energy, 
change in internal energy, work etc., in systems far from 
linear response regime into the relam of nonequilibrium 
domain. Some recent studies have also tried to reveal 
the relations between two completely unrelated phenom- 
ena, namely, stochastic resonance (SR) and Brownian 
ratchets in a formal way through the consideration of 
Fokker-Planck equations @, 0, [nl Q • It has been 
argued that the rate of flow of particles in a Brownian 
ratchet is analogous to the rate of flow of information 
in the case of stochastic resonance. Qian et al. have 
investigated a simple flashing ratchet model for ratchet 
effect as well as SR and have pointed out that the consis- 
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tency between these two phenomena are just due to the 
existence of circular flux in nonequilibrium state . 

In our present work we have analysed the nature of 
input energy (energy loss) and the total entropy pro- 
duction in a class of adiabatically rocked ratchets as a 
function of temperature of the bath (or noise strength). 
These systems exhibit peak in the noise induced current 
(in the absence of any net bias) as a function of tem- 
perature. The question now arises as to whether this 
peak is related to the underlying resonance due to the 
synchronization of the position of the particle with the 
external drive induced by the noise. Our analysis of in- 
put energy Ei n , rules out the presence of any resonance 
features in the dynamics of the position of the particle 
in these systems in the adiabatic regime. 

The presence of net currents in the ratchets increases 
the amount of known information about the system than 
otherwise. This extra bit of information comes from the 
negentropy or the physical information supplied by the 
external nonequilibrium bath. The amount of informa- 
tion transferred by the nonequilibrium bath is quantified 
in terms of algorithmic complexity. It has been argued 
that the algorithmic complexity or Kolmogorov informa- 
tion entropy is maximum when the current is maximum 
[l^ . Since the currents are generated at the expense of 
entropy we naively expect the maxmima in current to be 
related to the maxmima in the overall entropy produc- 
tion as a function of noise strength. However, we show 
that the maxima in current and the entropy production 
do not correlate with each other. 



II. THE MODEL: 

We study the motion of an overdamped Brownian 
particle in a potential V(q) subjected to a space de- 
pendent medium with friction coefficient j(q) and an 
external periodic force field F(t) at temperature T. The 
motion is described by the Langevin equation 



V'(q) - F(t) k B Ti(q) 



7(9) 



2[7(<?)r 



'k B T 
7(9) 



m (i) 



where £(t) is a randomly fluctuating Gaussian white 
noise with zero mean and correlation: < £(t)£(t') > 
= 2 S(t — t'). We take the potential V(q) to be periodic 
in space and is given by V(q) — ~sin(q) — (/i/4) sin(2q) 



with the asymmetry parameter u taking values between 
— I and 1. Also, we take the friction coefficient 7(g) to 
be periodic: 7(g) = 70(1 — Xsin(q + </>)), where (j) is the 
phase difference with respect to V(q) and the coefficient 
A takes values between and 1. The equation of motion 
is equivalently given by the Fokker-Planck equation 



dP{q,t) d I 



dt 



dqj(q) 



k B T 



dP(g,t) 
dq 



(2) 



[V'(q)-F(t)]P(q,t) 



This equation can be solved for the probability cur 
rent j when F(t) = Fq = constant, and is given by 

I - exp\- 27rFo] 

j 



k B T 



Jo* dyl-iy) 



(3) 



where 7_ (y) is given by 

i-{y) = 



exp 



-V(y) + F Q y 



fy +2 * dx 'f(x)exp 



k B T 
V(x) 



F x 



knT 



(4) 



In the case of inhomogeneous ratchets (space depen- 
dent frictional case, A / 0) flii \v\ with the spatial 
asymmetry parameter u = it may be noted that j(Fo) 
may not be equal to —j(—F ) for <f> ^ 0, tt. This fact 
leads to the rectification of current in the presence of an 
applied ac field F(t). In these inhomogeneous systems 
directed currents can be obtained even in a spatially pe- 
riodic symmetric potential. The inversion symmetry in 
these systems being broken dynamically by the space de- 
pendent frictional coefficient which is periodic in space 
having a phase lag of (f> with the potential profile. In 
the second case where A = (purely homogeneous case) 
the net currents are generated due to the spatial asym- 
metry of the potential (a ^ 0). We assume that F(t) 
changes slowly enough (adiabatic regime), i.e., its fre- 
quency is smaller than any other frequency related to 
the relaxation rate in the problem such that the sys- 
tem is in a steady state at each instant of time. For 
a field F(t) of a square wave of amplitude F , an av- 
erage current over the period of oscillation is given by, 
< j >= |L?(-Pb) +j(-F )} UlEi- In the quasi static 
limit following the method of stochastic energetics it can 
be shown [18| that the input energy Ei n (per unit time) 
is given by E m = ±F [j(F ) - j(-F )]. 

We also consider another type of ratchet, namely, the 
time asymmetric ratchets |2fJ, |2l| where the driving 
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force has zero mean, < F(t) > = 0, but is asymmetric 
in time i.e., 

F(t) = i±f [n T < t < nr + i T (i _ e )) 5 (5) 
= -F , (nr + -r(l - e) < t < (n + l)r). 
The time averaged current in this case is given by 

<3>=^{h+h) (6) 

with 

h = (l-e)j(^i^o) (7) 
h = (l + e)j(-F ) 

The input energy Ei n per unit time for this time asym- 
metric ratchet is given by E in = ^F a ( j^ji — 32) |21|. 
For the case where fi = 0, A = and e ^ the currents 
are generated in the absence of broken spatial symme- 
try, but in the presence of a temporal asymmetric driv- 
ing with zero mean. This type of temporal asymmetry 
is particularly common in biological systems [19j. 

III. RESULTS AND DISCUSSIONS 
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FIG. 1: El, < j > and S p vs temperature for n — 1.0 , A = 
0.0 , e = with fixed F = 0.5 and <j> = 0.3tt. 

In the following we analyse all these special classes of 
adiabatic ratchets mentioned above and also their com- 
binations. We study the average current < j >, the 
total entropy production S p , and the energy loss E L as 




Temperature 



FIG. 2: El, < j > and S p vs temperature for fj, = 0, A = 
0.9 , e = with fixed F = 0.5 and <j> = 1.3tt. 

a function of temperature T (noise strength) for seven 
different cases of adiabatically rocked ratchet systems 
described below. All these quantities are averaged over 
the period of external drive and are in appropriate di- 
mensionless units. It has been argued that for the case 
of a driven double well system input energy is a reli- 
able quantity for the identification of SR taking into 
account the detailed comparison between various mea- 
sures of SR |22j. Further analysis based on input en- 
ergy has shown the SR to be a bonafide resonance ji^f 
in that one obtaines peak in the input energy as a func- 
tion of noise strength as well as the frequency of the 
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Temperature 

FIG. 3: El, < j > and S p vs temperature for (j, = 0, A = 
, e = 0.6 with fixed F = 0.1 and (p = 0.3tt. 
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external drive Thus the peak in the input energy 
represents the matching condition of the escape rate out 
of the potential well (or synchonization in the dynam- 
ics of the particle) and the external driving frequency. 
It should be noted that as the system on the average 
does not perform any useful work the input energy in 
the steady state equals the energy loss (hysterisis loss) 
in the medium. Hysterisis loss being a good measure 
to identify SR is already known in the literature [2o|| . 
It is quite natural that when the system dynamics ex- 
hibits a resonance feature (or a peak) by tuning certain 
physical parameters, then at resonance the input energy 




Temperature 

FIG. 4: El, < j > and S p vs temperature for /i = 1.0 , A = 
0.9 , e = with fixed Fa = 0.5 and <p = 0.3vr 
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Temperature 

FIG. 5: El, < j > and S p vs temperature for /j, = , A = 
0.9 , e = 0.4 with fixed F = 0.3 and cj> = 0.3vr. 



extracted from the source (and the concomittant energy 
loss in the medium) is expected to be high. Thus the 
study of input energy or energy loss is expected to reveal 
the resonances if any in the dynamics of the particle as 
a function of various physical parameters. 

In our present work the particle performs a motion in 
a periodic potential in the presence of an adiabatic drive. 
As there is no load applied to the system, the system 
does not perform any useful work or stores energy and 

0.10 I ■ 1 . 1 1 1 . 1 1 1 
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-0.02 




-0.06 1 ' ' ' ' ' 1 ' 1 ' 1 

1 2 3 4 5 

Temperature 

FIG. 6: El, < j > and S p vs temperature for fi = —1.0 , A = 
0,e = 0.25 with fixed F = 0.6 and = 0.3tt. The current 
is scaled by a factor of 10 to make it comparable with El 
and S p . 




Temperature 

FIG. 7: El, < j > and S p vs temperature for n = —1.0 , A = 
0.9, e = 0.34 with fixed F = 0.3 and <j> = 1.005tt. The 
current is scaled by a factor of 10 to make it comparable 
with El and S p . 
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as a result all the input energy in the steady state will be 
dissipated away. Hence the energy loss in the medium 
is given by El, = Ei n . El inturn is equal to the heat 
dQ transferred to the bath and thus entropy production 
S p = dQ/T = E L IT 8]. Thus the total increase in the 
entropy (or the entropy production) of the bath (uni- 
verse) integrated over the period of the external drive 
is given by S p = ^ = ^ = ^ 0- As discussed 
in the introduction, currents in the ratchet systems are 
generated at the expense of entropy and thus we expect 
a correlation between the magnitude of current and the 
total entropy production. 

In figures 1 to 7 we have plotted energy loss El (equal 
to input energy), average current < j > and entropy 
production S p as a function of temperature T for var- 
ious types of adiabatically driven ratchet systems. All 
the physical parameters are in dimensionless units and 
their values are mentioned in the figure caption. All 
these figures are representative of the different classes 
of ratchet systems chosen to make our analysis clear. 

Figures 1,2 and 3 correspond to the spatially asym- 
metric case fi ^ 0, A = 0, e = 0; inhomogeneous case 
A 0, /i = 0, e = and temporal asymmetric case 
e ^ 0, A = 0, /i = respectively. In each case cur- 
rents are generated via different types of mechanisms 
arising out of spatial asymmetry in the potential or fric- 
tional inhomogenity or temporal asymmetry of the ex- 
ternal force. Here nonlinearity of the system, exter- 
nal drive and asymmetry conspire to generate unidi- 
rectional currents in the absence of any bias. In all 
these figures (1 to 3) current exhibits a single peak as 
a function of temperature. The input energy being a 
monotonically increasing function of T rules out any 
resonance in the dynamics of the particle as a func- 
tion of noise strength as discussed earlier. Entropy pro- 
duction exhibits a single peak as a function of noise 
strength. Moreover, the peak in the average current, 
< j > and total entropy production S p does not occur 
at the same T. This clearly indicates that maxima in 
the entropy production does not take place at the same 
value where current is maximum. To make this point 
explicit in figures 4,5 and 6 we have plotted Ei n , < j > 
and S p as a function of T. They correspond to a combi- 
nation of spatial asymmetry and system inhomogenity 
/i 7^ 0,A 7^ 0, e = 0; system inhomogenity and tem- 
poral inhomogenity A ^ 0, e ^ 0, /i = and temporal 



asymmetry and spatial asymmetry e^0,^^0,A = 
respectively. All these figures from 4 to 6 exhibit the 
phenomena of single current reversal or the absolute 
magnitude of current exhibits two peaks. However, the 
entropy production exhibits only a single maxima as a 
function of noise strength. This rules out clearly any 
correlation between current maxima and the maxima in 
the entropy production. The behaviour of input energy 
rules out any resonance phenomena in these systems as 
well. 

In fig. 7 we plotted < j > , Ej, n and S p for a sys- 
tem which incorporates frictional inhomogenity, tempo- 
ral asymmetry as well as spatial asymmetry. The sys- 
tem with all these combinations exhibit rich variety in 
the nature of current as a function of various physi- 
cal parameters and moreover in some parameter re gion 
large efficiency of energy transduction is observed 26]. 
For parameters given in fig. 7 the current exhibits two 
reversals as a function of T. Our observation of dou- 
ble current reversal in the adiabatic regime is in itself 
a novel phenomena. The absolute magnitude of cur- 
rent exhibits three peaks as against to the single peak 
structure in entropy production. This again reinforces 
the fact that maxima in the entropy production and 
current are totally unrelated. The behaviour of input 
energy again rules out the resonance dynamics in the 
system. 

In conclusion, by considering different cases of adia- 
batically rocked ratchets we have shown that the res- 
onance like feature observed in the nature of current 
as a function of temperature or the noise strength is 
not related to the intrinsic resonance in the dynamics of 
the particle and that the total entropy production does 
not extremize at the same parameter value at which the 
current exhibits a maximum. Our present results are 
valid only for the case of an adiabatically rocked ther- 
mal ratchet. This does not rule out the resonance in 
the dynamics of the particle in the nonadiabatic regimes 
as well as in other ratchet systems like flashing ratch- 
ets considered earlier |3, llfj, [llj, [l2( . The fact that the 
noise strength at which the current maxima and the 
maxima in entropy production occurs do not coincide 
may be related to the quality of the current. Noise in- 
duced currents are always associated with a dispersion 
or diffusion. When the diffusion is large then the quality 
of transport degrades and the coherence in the unidirec- 
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tional motion is lost. The coherent transport (optimal 
transport) refers to the case of large mean velocity at 
fairly small diffusion. It can be quantified 0,0001 
by a dimensionless Peclet number. For a given mag- 
nitude of current transport may be coherent or inco- 
herent. Thus analysis of the relation between current 
and the entropy production requires not only the mag- 
nitude of current but also the quality of transport. It 
is also shown that the algorithmic complexity or the 



Kolmogorov information entropy of the thermal ratchet 
motion exhibits maxima at the same value at which cur- 
rent is maximum |l4j | . It will be of interest to see mul- 
tiple maxima in algorithmic complexity as a function of 
system parameter in a multiple current reversal regime. 
These studies are expected to reveal a deep connection 
between efficiency, quality of transport, entropy and in- 
formation. Works along these different directions are in 
progress. 
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